Giant Resonances using Correlated Realistic Interactions: The Case for Second RPA 
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Lately we have been tackling the problem of describing nuclear collective excitations starting from 
correlated realistic nucleon-nucleon (NN) interactions. The latter are constructed within the Unitary 
Correlation Operator Method (UCOM), which explicitly considers short-range correlations in order 
to properly soften the short-range behaviour of realistic NN potentials. It has been concluded that 
first-order RPA with a two-body UCOM interaction (UCOM-RPA) is not capable, in general, of 
reproducing quantitatively the properties of giant resonances (GRs), due to missing higher-order 
configurations and long-range correlations as well as neglected three-body terms in the Ifamiltonian. 

In the present paper we report results on GRs obtained by employing a UCOM interaction, based 
on the Argonne V18 potential, in Second RPA (SRPA). The same interaction is used to describe 
the Hartree-Fock (HP) ground state and the residual interactions. We find that the inclusion of 
second-order configurations - which effectively dress the underlying HF single-particle states with 
self-energy insertions - produces sizable corrections. The effect appears essential for a realistic 
description of GRs when using the UCOM. We argue that effects of higher than second order should 
be negligible. Therefore, the UCOM-SRPA emerges as a promising tool for consistent calculations 
of collective states in closed-shell nuclei. This is an interesting development, given that SRPA can 
accommodate more physics than RPA (e.g., fragmentation). Remaining discrepancies due to the 
missing three-body terms and self-consistency issues of the present SRPA model are pointed out. 
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I. INTRODUCTION 



Many-body approximations like Hartree-Fock (HF) 
and Random-Phase Approximation (RPA) (and 
their counterparts for open-sheU nuclei, namely HF- 
Bogohubov (HFB) and Quasiparticle-RPA (QRPA)) 
have allowed massive calculations of nuclear ground- 
state and excited-state properties throughout the nuclear 
chart. Such models are used in conjunction with effective 
nucleon-nucleon (NN) interactions. Indeed, the bare NN 
interaction induces strong correlations in the nuclear 
system - most notably short-range correlations - which 
cannot be described by simple model spaces such as 
those involved in HF, RPA, etc. One can find very good 
parameterizations of the effective NN force - and there is 
now intense and coordinated activity towards the devel- 
opment of high-quality energy functionals to serve such 
a purpose - but those have been phenomenological up 
to now, lacking a direct connection with the underlying 
bare interaction. 

The question remains as to whether it is possible to 
construct a global effective NN interaction starting from 
the bare one. There have been two recent attempts to- 
wards that direction. One is the construction of a low- 
momentum interaction, the so-called T4ow-k, by integrat- 
ing out the high-momentum components of the bare one 
(thus softening its short-range behavior) using renormal- 
ization group techniques [l|. The other one is the Uni- 
tary Correlation Operator Method (UCOM) 0, S, Q, 
which deals explicitly with the short-range correlations 
and is described in the next Section. Applied to a realis- 



tic NN interaction, the UCOM produces a "correlated" 
interaction, Vucom- Although constructed following dif- 
ferent formalisms, 14ow-k and Vucom have similar low- 
momentum matrix elements.^ Moreover, they do not de- 
pend strongly on the particular bare NN potential on 
which they are based, and to which they are phase-shift 
equivalent. The hope is to be able to employ such realis- 
tic but "softened" potentials in many-body calculations. 

In the present work we focus on nuclear giant reso- 
nances (GRs) of closed-shell nuclei. First-order RPA with 
a two-body UCOM interaction has not been able to re- 
produce quantitatively the properties of all GRs 0, 0]. 
Here we report results on GRs obtained using Second 
RPA (SRPA) and employing the correlated Argonne V18 
interaction (UCOM-SRPA). 

It is not straightforward to perform SRPA calculations 
self-consistently - in the sense that exactly the same in- 
teraction is used to describe the ground state and the 
residual couplings - without conceptual problems. In 
typical SRPA applications in the past, phenomenologi- 
cal single-particle energies have been used and G-matrix 
or phenomenological forces have been employed as resid- 
ual interactions. (The real part of the SRPA self energy 
would then be discarded, since it would shift the already 
realistic single-particle energies.) Phenomenological den- 
sity functionals, on the other hand, are typically fitted by 
using HF(B) and (Q)RPA results. Part of the long-range 
correlations affecting ground-state properties are then ef- 
fectively taken into account by the parameterization and 
higher-order effects are usually ignored. Employing such 
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^ There remain important differences, however, as demonstrated, 
e.g., by the fact that the VijcOM is able to produce stable (sat- 
urated) nuclear matter, while Viow-k is not. 
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interactions in SRPA might result in overcounting of such 
effects.^ Our correlated interaction, however, takes into 
account only short-range correlations; long-range correla- 
tions have to be described by extending the configuration 
space and one way to do that is SRPA. 

Let us note at this point that the RPA reaches its 
limits when confronted with problems such as the width 
and fine structure of GRs, the strength of some low-lying 
states, etc. The SRPA is a more appropriate theory to 
deal with such issues. 

In the next section we outline the basic principles of 
the UCOM scheme. In Sec. IIIII we review what we have 
learned so far by using the Vucom in HF, perturbation 
theory, and first-order RPA calculations. In Sec. IIVI we 
present the SRPA formalism and our new results. In 
Sec. |V]we give a summary and perspectives. 



II. THE UCOM HAMILTONIAN 

The basic idea of the UCOM is the explicit treatment 
of the interaction-induced short-range central and ten- 
sor correlations. These are imprinted into an uncorre- 
lated many-body state \'^) (e.g., a Slater determinant) 
through a state-independent unitary transformation de- 
fined by the unitary correlation operator C, resulting in a 
correlated state \^) = C\'^). The correlation operator C 
is written as a product of unitary operators Cq and C,- 
describing tensor and central correlations, respectively. 
Both are formulated as exponentials of a Hermitian gen- 
erator, 

C = CfiCr = exp[-i^go,jj] exp[-i^gr,ij]. (1) 

The construction of the two-body generators gr and ga 
follows the physical mechanisms by which the interac- 
tion induces central and tensor correlations. The short- 
range central correlations, caused by the repulsive core 
of the interaction, are introduced by a radial distance- 
dependent shift pushing nucleons apart from each other 
if they are within the range of the core. Tensor corre- 
lations between two nucleons are generated by a spatial 
shift perpendicular to the radial direction. For a given 
bare potential, the corresponding correlation functions 
are determined by an energy minimization in the two- 
body system for each (S*, T) channel. 

Matrix elements of an operator O with correlated 
many-body states \^) can be equivalently written as 
matrix elements of a "correlated" (transformed) oper- 
ator O = C^OC and uncorrelated many-body states 



^ One could, in principle, consider to fit their parameters using 
SRPA, but that would be a formidable task from a computational 
point of view. Note that computationally friendly zero-range 
interactions are not appropriate for large-scale SRPA. 



1 4*). Thus, one can work in simple Hilbert spaces (sim- 
ple states) using correlated operators, rather than with 
bare operators and explicitly correlated states. By ap- 
plying the transformation to a bare NN interaction, a 
phase-shift equivalent correlated interaction is obtained, 
is suitable for use in tractable model spaces [3) 0, HI . The 
same transformation can then be applied to any other op- 
erator under study, as is needed for a consistent UCOM 
treatment. 

In an ^-body system a correlated operator contains 
irreducible contributions to all particle numbers. The 
cluster expansion of a correlated operator reads 

d = C^OC^ Ol^l + d[2l + ■ • • + Ol^l , (2) 

where O^"! denotes the irreducible n-body contribution. 
In actual applications of the UCOM a two-body approx- 
imation is usually employed, i.e., three-body and higher- 
order terms of the expansion are neglected. Starting from 
the uncorrelated Hamiltonian H for the v4-body system, 
consisting of the kinetic energy operator T and a two- 
body potential V, the formalism of the UCOM is used 
to construct the correlated Hamiltonian in two-body ap- 
proximation 

^C2^y[l]^j,[2]_^y[2] (3) 

where the one-body contribution comes only from the 
uncorrelated kinetic energy tI^I = T. Two-body contri- 
butions arise from the correlated kinetic energy f and 
the correlated potential V^'^\ which together constitute 
the phase-shift equivalent correlated interaction Vucom- 

It has been verified that higher-order contributions due 
to short-range central correlations can be neglected in 
the description of nuclear structure properties [J- The 
tensor interaction, on the other hand, is long-ranged and 
thus generates long-range correlations in an isolated two- 
nucleon system. However, the long-range tensor correla- 
tions between two nucleons embedded in a many-nucleon 
system are suppressed by the presence of other nucle- 
ons, leading to a screening of the tensor correlations at 
large interparticle distances. In order to effectively de- 
scribe the screening effect and at the same time justify 
the two-body approximation, the range of the tensor cor- 
relation function — more precisely, the "correlation vol- 
ume" I^)^''^^ — is restricted during the parameteriza- 
tion procedure. Restricting the range of the tensor corre- 
lator has another important function, namely to ensure 
that only state-independent, short-range correlations are 
described by the UCOM. By varying the correlation vol- 
umes — the only parameters entering the formalism — 
a family of correlators and respective correlated interac- 
tions are obtained. 

The question is then how to optimize these parameters 
in order to best describe the screening effect and the sep- 
aration of the two types of correlations. As demonstrated 
in Ref. this can be done with the help of exact few- 
body calculations. In particular, the values can be cho- 
sen so as to best describe the binding energies of "^H and 



3 



He within the no-core sheh model (4^'°^ = 0.09 fm^ for 
the Argonne V18 potential). For such a choice of tensor 
correlator range the missing genuine three-nucleon inter- 
action and the omitted higher-order terms of the cluster 
expansion of the correlated Hamiltonian effectively can- 
cel each other. As was subsequently shown within many- 
body perturbation theory f§\ , and verified by RPA calcu- 
lations [9], this cancelation remains at work throughout 
the nuclear chart, as far as the binding energy is con- 
cerned (see also next Section).'^ 

In this work we will use the correlated Argonne V18 
potential with 4^'°^ = 0.09 fm"^. No tensor correlator is 
employed in the triplet-odd channel, where the tensor in- 
teraction is much weaker. We start from a Hamiltonian 
which consists of the intrinsic kinetic energy Tint and the 
VucOM interaction derived from the Argonne V18 poten- 
tial including the Coulomb potential, 

-ffint = T — Tcin + VuCOM — Tint + ^UCOM , (4) 

in two-body approximation. It is the two-body Hamil- 
tonian Hint that has been used in Hartree-Fock (HF), 
perturbation-theory, and RPA calculations in Refs. [g, 
|8| and that will be employed in this work too. In prac- 
tice, two-body matrix elements in a harmonic-oscillator 
basis are the input to such calculations. 



III. APPLICATIONS IN SPHERICAL NUCLEI: 
RECENT LESSONS 

Using the Vucom in HF calculations we obtained 
bound nuclei throughout the nuclear chart Q. The ten- 
sor correlations play an important role in this. Note, 
though, that using the UCOM we aim to treat explic- 
itly only the state-independent short-range correlations; 
long-range correlations should be described by the model 
space. This tells us already that the UCOM-based HF is 
not enough, since a Slater-determinant wavefunction is 
unable to describe correlations. It is found indeed that 
the binding energies are underestimated by about 4 MeV 
per nucleon. The charge radii are underestimated too. 
While the Fermi energy is correctly reproduced, the level 
spacing of the single-particle states is too small. 

Second-order perturbation theory constitutes a 
tractable extension to the "zero-order" description 
provided by HF and was employed in Ref. [1]. The 
very good description of nuclear binding energies 
achieved within perturbation theory for nuclei from 
''He to ^°^Pb shows that the cancellation between the 
omitted three-body terms of the cluster expansion 
and genuine three-body correlations and terms of the 



The VljcOM has a strong momentum dependence, even when it 
is based on a local potential like Argonne V18. That is why it 
can perform reasonably well without an additional thrcc-body 
term. 



interaction works throughout the nuclear chart as far 
as the binding energies are concerned. Charge radii 
are still underestimated within perturbation theory, 
suggesting that the above-mentioned cancellation does 
not work for all observables and that supplementing our 
two-body Hamiltonian with a three-body term to take 
account of missing effects may be necessary for realistic 
nuclear-structure calculations. Higher than second-order 
corrections are found to be small. 

The Vucom has also been employed in standard, self- 
consistent RPA calculations to study nuclear giant res- 
onances Q. The ground state was described by the 
uncorrelated HF state, as usual. The isoscalar (IS) gi- 
ant monopole resonance (GMR), the isovector (IV) giant 
dipole resonance (GDR), and the IS giant quadrupole res- 
onance (GQR) were examined. Highly collective states 
were obtained for various closed-shell nuclei ranging from 
^^O to ^°*Pb. A reasonable agreement with the exper- 
imental centroid energies of the IS GMR was achieved. 
By contrast, the energies of the IV GDR and the IS GQR 
were overestimated by several MeV. 

Obviously, the Vucom is not a traditional effective in- 
teraction. Partly because no long-range correlations are 
(effectively) included in the UCOM, the corresponding 
nucleon effective mass in nuclear matter obtained in a 
HF calculation is very low (around half the bare nucleon 
mass). This is confirmed by the HF results in finite nu- 
clei, in particular the small level density. It is also mani- 
fested by the above-mentioned RPA results on the GQR 
and GDR centroids. It follows then that, besides the 
possible important role of missing three-body terms in 
the Hamiltonian, another source of our failure to describe 
nuclear collective states quantitatively with UCOM-RPA 
can be residual long-range correlations. 

The standard RPA is based on the assumption that the 
true RPA ground state can be approximated by the HF 
ground state. It is not obvious that this assumption holds 
when the Vucom is used, given the large correction to the 
HF binding energies due to second-order ^ and RPA P 
correlations. Therefore, in Ref. @ the effect of explicit 
RPA ground-state correlations on the results for GRs was 
examined. To this end, a renormalized RPA version was 
used [13, UH El]- The effect on the properties of GRs 
was found to be rather small. It is concluded that first- 
order RPA with a two-body UCOM Hamiltonian cannot 
describe quantitativly the properties of GRs. 

Up to now we have assumed that residual three-body 
forces can be neglected, based upon the fact that they 
contribute only marginally to the ground-state energy as 
calculated within many-body perturbation theory 0, • 
This is not necessarily a valid assumption. A simple phe- 
nomenological zero-range three-body force can be con- 
structed in order to be used along with the correlated 
two- nucleon interaction in future calculations. Prelimi- 
nary results show that by using such a three-body force 
it is possible to improve on the description of observables 
such as nuclear radii and resonance energies while retain- 
ing the good reproduction of the experimental binding 
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energies. 

Another important issue with RPA is that only one- 
particle-one-hole excitations are taken into account and 
the coupling to higher-order configurations {2p2h and be- 
yond) is neglected. One can include higher-order con- 
figurations, starting with two-particle-two-holc within 
SRPA. Given that an extended model space is of great 
importance when using the Vucom, it is imperative to 
examine the effect. 



IV. SECOND RPA 
A. Formalism 

We will use the SRPA as it was formulated in Ref. 
in analogy to RPA. Excited states Iv) of energy E^, = %uiv 
with respect to the ground state |0) 



IW = Qt|o), 



(5) 



are considered as combinations of Iplh and 2p2/i con- 
figurations. (We omit angular momentum coupling to 
keep the notation simple.) The corresponding creation 
operators are then written as 



~''Z-ipihip2h2 Plhip2h2^pihip2h2 



(6) 



where O^^ creates a ph state and 0^,^^,^, creates a 2p2h 
state. The SRPA ground state, which is the vacuum 
of the annihilation operators Q^, is approximated by the 
HF ground state. The forward {X, X) and backward {Y , 
y) amplitudes are then given by the SRPA equations in 
ph © 2p2 ft,— space 
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where A and B are the usual RPA matrices, A\2 describes 
the coupling between ph and 2p2h states and A22 con- 
tains the 2p2h states and their interactions. If we neglect 
the coupling amongst those states, A22 is diagonal and 
its elements are equal to the unperturbed 2p2h energies, 

A22 — Spip'^Shih{5p^p'^Shih[iepi + — e^^ — e/12) (8) 

{ci are the HF single-particle energies). 

The dimension N of the SRPA matrix, Eq. ([7]), can be 
very large. Fortunately, the SRPA matrix is also sparse. 
When the approximation ^ is employed, most of its 
elements are zero. Thus it is possible to store all its 
finite matrix elements and then use a Lanczos procedure 
to obtain only the eigenvectors of interest. 

The SRPA problem, Eq. ([7]), can be reduced to an 
energy-dependent eigenvalue problem of the dimension 



of the RPA matrix (see, e.g., Ref. [ij)- Therefore, it can 
be viewed as an RPA problem with an energy-dependent 
interaction. In general, the reduction procedure involves 
the inversion of a complex matrix in the large 2p2h space, 
but when A22 is diagonal, that is reduced to a trivial 
complex-number inversion. There are ways to solve an 
energy-dependent eigenvalue problem [Tsl . [iGj . An effi- 
cient alternative is to employ the response-function for- 
malism. Then, instead of explicitly solving the eigenvalue 
problem, one can obtain directly the strength distribu- 
tion of interest [3, . We have used this technique as 
well. 

It has been shown formally that the total strenght toq 
and the first moment of the strength distribution mi are 
the same in the present SRPA as in RPA fl7^. However, 
when based on the HF ground state, the SRPA is not 
fully self-consistent and symmetry-conserving, contrary 
to the RPA based on the HF ground state. It has been 
pointed out [18,] that it misses a class of second-order ef- 
fects, related to ground state correlations. The missing 
effects may be important, especially for the less collective 
low-lying states. In principle, it is possible to combine 
the SRPA with a correlated ground state [H, [l^, H^l for 
a most complete theoretical treatment of nuclear exci- 
tations, but that is beyond the purposes of the present 
work. 



B. Results 

We have used the correlated Argonne VIS interaction 
and a single-particle basis of 11 oscillator shells and we 
have examined the IS monopole (ISM), IV dipole (IVD) 
and IS quadrupole (ISQ) response of the nuclei ^^O and 
^"Ca. The convergence of the GR sum rules mo and mi 
and centroids is rather good for the present basis (within 
about 1 MeV for the centroids). The total number of 
eigenvalues is 10'''"^ for the cases presented here (it can 
be larger e.g. for heavier nuclei), but less than 300 eigen- 
states are sufficient to describe the region of the GRs. 
We use standard single-particle transition operators 
We present our results in comparison with experimental 
data. The experimental centroids mi /mo of the IS GMR 
and the IS GQR were taken from Refs. i2l| (i^O) and 
[2^ C^^Ca). Photoabsorption cross sections were found 
in Refs. Ij,!!! (^^O) and ^ {^°Ca.)^ and strength dis- 
tributions and centroids of the IV GDR were evaluated 
from those. 

We have verified that in SRPA the total mo is almost 
the same as in RPA. The total mi is smaller by more 
than 20%. This is probably because we do not calcu- 
late the full spectrum. A non-negligible part of the total 
mi may be distributed among the large number of weak 
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excitations lying at high energies. 

It has been shown [2^ that the spurious state related 
to the CM motion will generally not be exactly seper- 
atcd from the physical spectrum, when SRPA is based 
on the HF ground state. In order to quantify this prob- 
lem, we have examined the behaviour of the IS dipolc 
response. We found that the spurious state appears at 
about 5 MeV. We used a transition operator of the usual 
radial form (oc — ^{r'^)r) and its uncorrected form 
(oc r^) and found that the spectrum beyond the spuri- 
ous state is practically the same and can be considered 
uncontaminated. 

In Fig. [1] we show the ISM, IVD and ISQ strength 
distributions for the two nuclei. Note that, for presen- 
tation purposes, the calculated distributions (RPA and 
SRPA) have been folded with a Lorenzian with a width of 
2 MeV. Thus, all peaks have acquired an artificial width 
(which for some low-lying dipole states may be too large) 
and the SRPA fragmentation is not visible. In all cases, 
the SRPA centroid energies are much lower than the RPA 
ones. The reason for the difference between the RPA and 
SRPA results ~ even for such collective ph excitations like 
the GRs considered here - is to a large extent that, within 
SRPA, the coupling of single-particle states with virtual 
phonons is implicitly taken into account. The inclusion 
of second-order configurations within SRPA effectively 
dresses the underlying HF single-particle states with self- 
energy insertions and brings them closer to each other 
energetically, thereby lowering the underlying ph ener- 
gies. It is an important physical effect which cannot be 
ignored when using completely "undressed" (with respect 
to long-range correlations) HF states like the ones pro- 
duced by the Vucom- In this scheme the undressed HF 
energies are viewed as auxiliary model quantities which 
should not be directly compared with experiment. 

Let us look at the results in more detail. In the middle 
panels of Fig.[T]the RPA and SRPA strength distributions 
are shown for the IV GDR, along with those extracted 
from experimental data (there has been no ad hoc renor- 
malization imposed). We observe that the IV GDR is 
more realistically reproduced within SRPA than within 
RPA. Its centroid energy is somewhat underestimated. 
In the lower panels we show the ISQ strength distribu- 
tions. The RPA and SRPA results are shown and the 
experimental centroids of the IS GQR are indicated. The 
agreement of the SRPA results with experiment is very 
good. It appears as though, once coupling to higher-order 
configurations is taken into account, a realistic effective 
mass is restored. 

In the upper panels of Fig. [1] we show the ISM strength 
distributions. The energies of the IS GMR are underesti- 
mated within SRPA. This is another indication that there 



are missing three-body effects and our two-body interac- 
tion should be supplemented with a three-body term to 
describe them. Normally, residual three-body corrections 
should affect the IS GMR most of all, since it is a com- 
pression mode. They should affect less strongly the IV 
GDR, where the nuclear interior plays a lesser role, and 
less the IS GQR, which is a surface mode. These physical 
arguments could serve as a guide for the construction of 
an appropriate effective three-body term. 

In the above, the approximation ([5]) has been used. It 
has been verified that inclusion of the couplings amongst 
the 2p2h states produces negligible corrections. Note 
that those couplings constitute higher-order effects. The 
indications that we had from our perturbation-theory re- 
sults, that corrections beyond second order are small, are 
thus confirmed. 



V. SUMMARY AND PERSPECTIVES 

We have employed a correlated interaction Vucom 
based on the Argonne V18 potential in SRPA calcula- 
tions of nuclear GRs. Short-range correlations are ex- 
plicitly taken into account. The same interaction is used 
to describe the Hartree-Fock (HF) ground state and the 
residual interactions. We found that the second-order 
configurations produce sizable corrections with respect 
to first-order RPA. They do so by effectively dressing 
the underlying HF single-particle states with self-energy 
insertions. The effect appears essential for a realistic 
description of GRs when using the Vucom- Effects of 
higher than second order should be negligible. Therefore, 
the UCOM-SRPA model emerges as a promising tool for 
consistent calculations of collective states in closed-shell 
nuclei. This development is interesting, given that SRPA 
can accommodate more physics than RPA (e.g., fragmen- 
tation width and fine structure of GRs). Remaining dis- 
crepancies, regarding in particular the IS GMR, can be 
attributed to missing three-body effects. Self-consistency 
issues of the present SRPA formulation were also pointed 
out. 

Up to now we have considered mostly the centroids 
of GRs, but their decay properties can also be studied 
within UCOM-SRPA. Heavier nuclei and low-lying states 
will be a topic for future work as well. 
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